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ABSTRACT 

We investigate the synchrotron self-Compton process in a planar shell taking 
the shock structure into account. We find that the energy density of the seed 
photons could deviate from the one-zone estimate by order of unity depending 
on the shock velocity and the electron cooling time. We also find that as the 
electron cooling becomes faster, the seed photons are increased more, so that the 
inverse Compton cooling becomes more efficient. This "ultra" fast cooling may 
work in such as gamma-ray bursts, blazars and microquasars. 

Subject headings: galaxies: jets — gamma rays: bursts — gamma rays: theory 
— radiation mechanisms: non-thermal — shock waves 



1. INTRODUCTION 

Relativistic shocks often arise in astrophysics when a faster flow hits upon a slower 
one such as in gamma-ray bursts (GRBs) (e.g., Piran 1999), blazars (e.g., Inoue & Taka- 
hara 1996; Kino, Takahara & Kusunose 2002; Rees 1978), microquasars (e.g., Mirabel & 
Rodriguez 1999; Levinson & Waxman 2001; Kaiser, Sunyaev & Spruit 2000) and so on. In 
the relativistic shocks, the kinetic energy of the flow turns into the internal one, and some 
fraction of the internal energy is distributed to electrons and magnetic fields. The electrons 
are accelerated in the shock front while the magnetic fields are amplified by the shock. Under 
these conditions, the accelerated electrons radiate nonthermal emission, such as synchrotron 
and inverse Compton (IC) emission. In this paper we will consider the synchrotron self- 
Compton (SSC) process, in which the seed photons for the IC emission are the synchrotron 
photons. 

The intensity of the IC emission is proportional to the energy density of the seed photons 
(Rybicki & Lightman 1979). So far the energy density of the seed photons has been estimated 
from the one-zone argument (e.g., Inoue & Takahara 1996; Sari, Narayan & Piran 1996; Sari 
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& Esin 2001). In the one-zone model we assume that the seed photons are distributed 
uniformly over space. However the one-zone approximation is too crude at times. In fact the 
shock structure, such as the cooling layer (e.g., Granot, Piran & Sari 2000) or the magnetized 
layer (e.g., Rossi & Rees 2002), could bring new features. Also for the seed photons, the 
shock structure may be important since the uniformity is violated near the shock front. At 
least we should evaluate the validity of the one-zone approximation quantitatively taking 
the shock structure into account. 

In this paper, we will investigate the effects of the shock structure on the energy density 
of the seed photons in the SSC process. Especially we will concentrate on the geometrical 
effects. We will find that for some parameters the energy density of the seed photons could 
deviate from the one-zone estimate by order of unity. 



Let us consider the following simple model. Our model is not so sophisticated as the 
realistic ones. However simplicity should be helpful in elucidating the main features. 

We consider an optically thin uniform shell with a thickness D and a shock propagating 
with a velocity c(3 = c(l — 7~ 2 ) 1//2 measured in the shocked fluid frame (see Fig. 1). Electrons 
are accelerated just behind the shock and the motion of the electrons is negligible in the 
shocked fluid frame, if the electron acceleration time is shorter than other timescales. To 
extract the geometrical effects, we assume that all particles and magnetic fields are isotropic 
in the shocked fluid frame. We consider only single scattering assuming that the higher order 
IC is suppressed by the Klein- Nishina effect. 

Note that if the shock is strong and the shocked fluid is extremely hot, the shock velocity 
measured in the shocked fluid frame is given by (3 = [(7,. — l)/(7r + l)] 1 ^ 2 /3 (Blandford & 
Mckee 1976), where 7 r ~ {if /is + 7s/7/)/2 is the relative Lorentz factor between the faster 
fluid with the Lorentz factor 7j and the slower fluid with j s . When 7//7 s ~ 2 and 7//7 s 1, 
P ~ 0.1 and (3 = 1/3, respectively. Therefore we will consider the cases (3 = 0.1 and (3 = 1/3. 

The energy density of the seed photons at a point r = R is given by the integration of 
all synchrotron radiation from the accelerated electrons as 



where P{t,r) erg s x cm 3 is the synchrotron power at a time t and a position r, \x = cos 8 
and s > (see Fig. 1). Note that the retarded time t — s/c is used in equation (1). For 
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PLANAR MODEL 
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simplicity, we consider the monoenergetic injection of the electrons. Since we concentrate on 
the total energy density of the seed photons, we can regard these monoenergetic electrons as 
the electrons that contribute most of the radiation energy even for the power-law injection 
case. To extract the essence, we adopt the following simple form as the synchrotron power, 



P(t,r) = A 



H [t — 



c/3 



— H [t — At coo i 
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(2) 



which means that electrons radiate with a constant power A erg s _1 cm~ 3 for a duration 
Atcooi s since the shock has passed. This is sufficient for the following argument, since the 
synchrotron power is proportional to the square of the electron Lorentz factor (Rybicki & 
Lightman 1979) and hence the cooled electrons have a small contribution to the total energy 
density of the seed photons. 

Equation (2) makes it possible to integrate equation (1) analytically. As shown in § A, 
during the electron emission, 
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the energy density in equation (1) is given by 
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are satisfied. If equations (3) and (6) are satisfied but equation (5) is not, we have 



U,(t,R) = ^ 
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(3) 



(4) 



(5) 



(6) 



(7) 



3. ENERGY DENSITY OF SEED PHOTONS 

Let us consider the energy density of the seed photons when the shock is in the region 
< r < D. The shock crossing time is given by At dyn = D/c(3. We call the case / = 
At dyn / Atcooi > 1 fast cooling, and f < I slow cooing. 
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First, let us consider the one- zone model. In fast cooing, electrons emit almost all energy 
before the shock crosses the shell, so that the energy density of the seed photons is given by 



lzone 
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dtP(t,r) = AAt cool , 



(8) 



with equation (2). In slow cooing, electrons do not emit all the energy within Atd yn - The 
electrons radiate with a power A for a duration Atdyn, so that 
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In the planar model, the energy density of the seed photons is given by equations 
(4) and (7). To compare the planar model with the one-zone model, we take the time 



and position average of the energy density, {U 7 (t,r)) = J Q (dr/D) f (dt/ At)U 7 (t,r) where 
At = min[At coo ;, (D — r)/cf3\. In the fast cooling limit / 3> 1, we have 



(U y (t,r)) = 
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Note that the logarithmic term originates from the integral of ~ i.e., the geometrical 
effect (see § A). In the slow cooling limit / < 1, we have 



(U y (t,r)) = 
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In Fig. 2 we show the ratio of the energy density of the seed photons to the one- 
zone estimate (U 7 (t, r))/U} / zone for (3 = 1/3 and f3 = 0.1. The analytical approximations 
in equations (10) and (11) are also shown by dashed lines. We see that the analytical 
approximations are quite good. 

From Fig. 2, we can find the following features. First, for some parameters, e.g., f3 = 0.1 
and / < 1, the energy density of the seed photons deviates from the one- zone estimate by 
order of unity. Second, the dependence of the energy density on the shock velocity cf3 differs 
between in the planar model and the one-zone model. Finally, in fast cooling the energy 
density of the seed photons increases as the ratio / increases, because of the logarithmic 
term in equation (10). This is an interesting new feature due to the geometrical effect. 
Increasing the seed photons enhances the IC emission and hence the IC cooling. Therefore 
this mechanism may be termed "ultra" fast cooling. The difference between / = 1 and 
/ ~ 10 6 is about order of unity when AAt coo i = const. Surprisingly the energy density of 
the seed photons can be larger than that of the electrons which supply the seed photons. 
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To make the physical situation clear, the region from which photons come is shade in 
Fig. 3. Here we assume that equations (3), (5) and (6) are satisfied. This causal region is 
bounded by four constraints as in Fig. 4 (see § A), (a) s = R/[i, (b) s = —(D — R)/fi, (c) 
s = (cf3t — R) I {j3 — /i) and (d) s— [cj3{t — At cool ) — R]/({3 — /i). The boundary (c) represents 
the hyperboloid of two sheets, 

x 2 + y 2 (r-R-^cPt) 2 _ 

7 2 c 2 f3H 2 7 4 c 2 /3 4 t 2 ' 1 ' 

where x 2 + y 2 = s 2 (l — /i 2 ). The boundary (d) is also described by equation (12) with t 
replaced by t — At coo i. From equation (12) we can see that the photons approximately come 
from the direction ~ tan 6* = / 9~ 1 7~ 1 . This is physically reasonable since the cooling layer 
for fast cooling is very close to the shock front and the shock front has the same velocity 
towards the r-axis as the photons traveling in the direction of tan# = / 3~ 1 7 _1 . In a sense, 
photons are accumulated around the shock front as the shock sweeps. The approximation 
of the planar shell is valid when the size of the shell is larger than ~ j3~ 1, -f~ 1 D. 



4. APPLICATIONS 

(I) In the internal shocks of the GRBs (e.g., Piran 1999; Sari, Narayan & Piran 1996), 
the cooling time of electrons is about At coo i ~ QTi^ e m e c 2 /cr T c-f 2 B 2 ~ 10~ 6 7 e " 3 1 J B 2 s, where 
7 e = 10 3 7 e 3 ~ m p /m e is the Lorentz factor of the electrons which contribute most of the 
radiation energy, and B = 10 6 i?6 G is the magnetic field. The shock crossing time is about 
Atdyn ~ Td/c ~ lT 2 d 8 s, where T = 10 2 T2 is the Lorentz factor of the shell and d = 10 8 d 8 cm 
is the shell thickness in the lab frame. A large dispersion of the Lorentz factors may be needed 
for the efficient internal shock (Beloborodov 2000; Kobayashi & Sari 2001). Thus f3 ~ 1/3 
and / = At dyn / At coo i ~ 10 6 . From Fig. 2, the energy density of the seed photons is larger 
than the one-zone estimate by a factor of ~ 3. 

(II) The late afterglows of the GRBs are likely in the slow cooling regime (Sari & Esin 
2001). Most of the radiation energy is emitted by electrons with / = 1. For the forward 
shock the shock velocity is (3 ~ 1/3, while for the reverse shock the shock velocity has various 
values depending on the parameters (Sari & Piran 1995). From Fig. 2, the one-zone estimate 
may overestimate the energy density of the seed photons by order of unity for (5 — 0.1. 

(III) In the TeV blazars, most of the radiation is emitted by electrons with the maximum 
Lorentz factor ^y max ~ 10 5 if the power law index of the electron distribution is s < 2 (see 
Kino, Takahara & Kusunose 2002). Since electrons with the break Lorentz factor / ~ff )r ~ 10 3 
corresponds to the electrons with / = 1 and the cooling time is proportional to the inverse 
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of the electron Lorentz factor, the electrons with 7 max have / ~ 10 2 . If we take the Klein- 
Nishina effect into account, photons emitted by the electrons with , ~f rnax cannot be seed 
photons for the electrons with 7 max . In this case / is reduced to / ~ 10 or so. If the 
emission arises from the internal shock with the relative Lorentz factor 7 r ~ 2 (Rees 1978), 
j3 ~ 0.1. From Fig. 2, the one-zone estimate may overestimate the energy density of the seed 
photons by a factor of ~ 4. 

(IV) In the microquasars, according to Levinson & Waxman (2001), the ratio of the 
shock crossing time to the cooling time is about / ~ lO 5 /^ 1 where / is the collision radius of 
the internal shock. If we take the Klein-Nishina effect into account in the case Is — 1, / is 
reduced to / ~ 10 4 or so. If the relative Lorentz factor between shells is 7 r ~ 2, (3 ~ 0.1. 
From Fig. 2, if / = 10 8 cm, the one-zone estimate is correct within a factor of 2, while if 
I = 10 13 cm, the one-zone estimate may overestimate the energy density of the seed photons 
by order of unity. 

5. DISCUSSIONS 

We have not dealt with the emission spectrum since we concern the total energy density 
in this paper. If we consider the energy distribution of electrons, the ultra fast cooling may 
modify the emission spectrum conventionally used in the one-zone model. This is because 
the cooling time depends on the electron energy, i.e., less energetic electrons cool slower. 
As expected from our analyses, the seed photon field that electrons feel will depend on the 
cooling time of the electrons. Therefore the spectrum of the IC emission should be modified, 
since the seed photon density is independent of the electron cooling time in the one-zone 
model. Furthermore if the IC cooling dominates the synchrotron one, the well known fact 
that the energy distribution of the cooled electrons breaks by one power in the index (e.g., 
Heavens & Meisenheimer 1987) should be modified because of the same reason. Detail 
studies will be presented in the future paper. 

To be precise, when the IC cooling dominates, we have to solve the electron cooling 
together with equation (1). This problem is the one-dimensional radiative transfer with the 
electron cooling. We may be able to solve this problem numerically, although a large memory 
will be necessary for resolving the cooling layer in fast cooling. 

We may have to consider other effects. For example, the approximation of the uniform 
shell may not be good. The reverse shock emission and the deceleration of the shock may 
be important. If the synchrotron power is not constant but a function of the distance from 
the center, the logarithmic correction may become a power law one. These are interesting 
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future problems. 
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In Fig. 4, the region to be integrated in equation (1) is shaded. There are four constraints 
on the integral in equation (1), which are shown by solid lines in Fig. 4. (a) s < Rj fj, since 
the shell has an end r = R — s/i > and we consider R > 0. (b) s < —(D — R)/ /j, since the 
shell has an end r — R — s/j, < D and we consider R < D. (c) s < {cf3t — R)/(/3 — /J,) since the 
first Heaviside step function in equation (2) has to be unity at the retarded time t — s/c, i.e., 
t — s/c—(R—sfi)/c(3 > 0. (d) s < [c/3(t — Atcooi) — R]/ (/3 — [/,) since the second step function in 
equation (2) has to be zero at the retarded time t — s/c, i.e., t — At coo i — s/c—(R — sn)/c(3 < 0. 
Here we consider the interval during which electrons are emitting in equation (3), so that 
the line (c) is in the region /i < (5 and the line (d) is in /i > f3 as in Fig. 4. The lines (a) and 
(d) cross in the region (3 < fi < 1 as in Fig. 4 if equation (5) is satisfied, and vice versa. The 
lines (b) and (c) do not cross in the region — 1 < /i < as in Fig. 4 if equation (6) is satisfied, 
and vice versa. If equations (3), (5) and (6) are satisfied, equation (1) can be integrated as 



which gives equation (4). If equations (3) and (6) are satisfied but equation (5) is not, we 
have 



A. INTEGRATION IN EQUATION (1) 
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which gives equation (7). 
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Fig. 1. — Our simple model is shown. We consider an optically thin uniform shell with a 
thickness D and a shock propagating with a velocity c(3 — c(l — 7 -2 ) 1 ^ 2 measured in the 
shocked fluid frame. Electrons are accelerated just behind the shock, and then they cool 
radiating synchrotron and Compton emission. The motion of the electrons is neglected in 
the shocked fluid frame. 
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Fig. 2. — The energy density of the seed photons divided by the one-zone estimate 
(£/ 7 )/£/,J zone as a function of the ratio of the shell crossing time to the electron cooling 
time / = Atdyn/ At coo i = D/c(3At coo i is shown for the shock velocities (3—1/3 and (5 — 0.1 
with solid lines. The analytical approximations in equations (10) and (11) are also shown 
by dashed lines. We call f > I fast cooling and / < 1 slow cooling. 
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Fig. 3. — The region from which photons come to the square point is shaded (see also Fig. 1), 
where we assume that equations (3), (5) and (6) are satisfied. This causal region is bounded 
by four constraints as in Fig. 4, (a) s = Rj fi, (b) s = —(D — R)/fi, (c) s — (c/3t — R)/(f3 — /i), 
and (d) s = [cf3(t — At coo i) — R]/(f3 — fi). The boundaries (c) and (d) are the hyperboloid of 
two sheets. The boundary (b) is not relevant to this figure. Approximately, photons come 
from the direction tan# = / 3~ 1 7~ 1 , which is shown by the dot-dashed lines. 
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Fig. 4. — The region to be integrated in equation (1) is shaded in the (//, s) plane. Here 
we assume that equations (3), (5) and (6) are satisfied. Four constraints on the integral are 
shown by solid lines, (a) s = Rj /i, (b) s = —(D — R)/ fi, (c) s = (cj3t — R)/(fl — /i) and (d) 
s = [c/3(t — At coo i) — R]/(f3 — fi). The lines (a) and (c) cross at /i = R/ct, and the lines (a) 
and (d) cross at ji = R/c(t — At coo i). 



